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Pulsars are known to power winds of relativistic particles that can produce 
bright nebulae by interacting with the surrounding medium. These pulsar 
wind nebulae (PWNe) are observed in the radio, optical, x-rays and, in some 
cases, also at TeV energies, but the lack of information in the gamma-ray band 
prevents from drawing a comprehensive multiwavelength picture of their phe- 
nomenology and emission mechanisms. Using data from the AGILE satellite, 
we detected the Vela pulsar wind nebula in the energy range from 100 MeV 
to 3 GeV. This result constrains the particle population responsible for the 
GeV emission, probing multivavelength PWN models, and establishes a class 
of gamma-ray emitters that could account for a fraction of the unidentified 
Galactic gamma-ray sources. 

The Vela supernova remnant (SNR) is the nearest SNR (d ~ 290 pc) containing a bright 
pulsar, PSR B0833-45, which has a characteristic age of 1 1 kyr, and a spin-down luminosity of 
7 x 10 36 erg s _1 ([7] |2j). This SNR extends over a diameter of ~8°, and is known from early radio 
observations to embrace a number of regions of non-thermal emission © including Vela X, a 
lOO'-diameter, flat-spectrum radio component near the center of the SNR. Vela X, separated by 
~40' from PSR B0833-45, is generally interpreted as the pulsar's radio synchrotron nebula (0 
[5]). A diffuse emission feature (~1° long) coincident with the centre of Vela X was detected 
in x-rays (0.6-7.0 keV) by the Rontgen © and ASCA © satellites. It was first suggested that 
this feature, which is closely aligned with a filament detected at radio wavelengths, corresponds 
to the outflow jet from the pulsar's pole ©. More recently, observations with Chandra © 
clearly unveiled the torus-like morphology of the compact x-ray nebula surrounding the pulsar 
and indicated that the centre of Vela X lies along the extension of the pulsar equator, although 
bending to the southwest. 

The detection of very high-energy (VHE; 0.5-70 TeV) gamma-rays from the Vela X region 
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was claimed by HESS (E/J and confirmed by CANGAROO CLD- The strong VHE-source, 
HESS J0835-455 (luminosity of ~10 33 erg s _1 at energies above 0.55 TeV), coincides with 
the region of hard x-ray emission seen by the Rontgen Satellite. The best-fit VHE emission 
centroid (Ra. = 08 h 35 m l s , Decl. = -45°34'40") is ~0.5° from the pulsar position and the 
VHE emission has an extension of ~5x4 parsec 2 . The detection of Vela X at TeV energies 
demonstrated that this source emits non-thermal radiation, in agreement with the hypothesis 
that it corresponds to the pulsar wind nebula, displaced to the south by the unequal pressure of 
the reverse shock from the SNR (1721) . 

The multiwavelength spectrum of the center of Vela X can be modeled as synchrotron ra- 
diation from energetic electrons within the cocoon (radio and x-rays) and inverse Compton 
emission from the scattering (by the same electron population) of the cosmic microwave back- 
ground radiation (CMBR), the Galactic far-infrared radiation (FIR) produced by reradiation of 
dust grains, and the local starlight (1731 U4\ [771) . Alternatively, a hadronic model can be invoked 
for the gamma-ray emission from the Vela X cocoon, where the emission is the result of the de- 
cay of neutral pions produced in proton-proton collisions (TioT) . Observations in the MeV-GeV 
band (HE) are crucial to distinguish between leptonic and hadronic models as well as to identify 
specific particle populations and spectra. 

The Vela region was recently observed from 30 MeV to 50 GeV by the AGILE (TTTI) and 
Fermi ([781) gamma-ray satellites. The Vela pulsar is the brightest persistent source of the GeV 
sky and, due to the limited angular resolution of the current-generation gamma-ray instruments, 
its gamma-ray pulsed emission dominates the surrounding region up to a radius of about 5°, 
preventing the effective identification of weaker nearby sources. 

The AGILE satellite (1791) observed the Vela pulsar for ~1 80 days (within 60° from the cen- 
ter of instrument's field of view) from 2007 July (54294.5 MJD) to 2009 September (55077.7 
MJD). To obtain precise radio ephemeris and model the Vela pulsar timing noise for the entire 
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AGILE data span, we made use of observations with the Mount Pleasant radio telescope. The 
Vela pulsar timing analysis provided a total of ~40000 pulsed counts with energies between 30 
MeV and 50 GeV; the difference between the radio and gamma-ray ephemeris was <lCT n s. 
Gamma-ray pulsed counts are concentrated within the phase interval 0.05-0.65 (where is the 
phase corresponding to the main radio peak). We verified that no pulsed gamma-ray emission 
is detected outside this interval, as reported by previous observations with EGRET (1201 1271) . 
AGILE (UJ) and Fermi (UJ). 

With the aim of performing a sensitive search for close faint sources excluding the bright 
emission from the Vela pulsar, we discarded the time intervals corresponding to the phase inter- 
val 0.05-0.65. The analysis of the resulting off-pulse images (taking only events corresponding 
to 0.65-1.05 pulsar phase interval, for a total of ~ 14000 events), unveiled few gamma-ray 
sources, none of which coincides with the Vela pulsar. A maximum likelihood analysis (1791) . 
performed on the E > 100 MeV dataset within a region of 5° around the pulsar position, re- 
vealed two sources at better than 3a confidence (see Figure 1): AGL J0848-4242 (at Galactic 
coordinates, I = 263.11°, b = 0.65°, 68% confidence error circle radius ~0.25°) and AGL 
J0834-4539 (at I = 263.88°, b = -3.17°, e.c. radius ~0.20°). A gamma-ray source coinci- 
dent with the EGRET source 3EG J084 1-4356 (1221) was also detected with lower significance, 
and the Vela Junior (RX J0852.0-4622) SNR (|25ll2?l) is also possibly contributing to a counts 
excess in the Galactic Plane around /~265.6°. 

The brightest gamma-ray source, AGL J0834-4539 (~5.9a significance, ~264 counts, 
F7 = (35 ± 7) x 10~ 8 ph cm" 2 s" 1 at E > 100 MeV), is located ~0.5° southwest from 
the Vela pulsar position (outside the 95% source position confidence contour) and has a spatial 
extent of ~1.5xl square degrees. Its shape is asymmetric and incompatible with the AGILE 
point-spread function. Therefore, possible residual emission from the pulsar (in principle asso- 
ciated to undetected weak peaks in the off-pulse interval of the light curve) cannot substantially 
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contribute to this diffuse feature. No relevant systematic errors on positions, fluxes and spectra 
(mostly due to uncertainties on the Galactic gamma-ray diffuse emission model) affect AG- 
ILE sources detected around 5a level. AGL 10834^4-539 is positionally concident with HESS 
J0835-455, the TeV source that is identified with the Vela X nebula, and has a similar brightness 
profile to it (Figure 1). This implies that AGL J0834-4539 is associated with the pulsar's PWN. 

Based on the available count statistics, we performed a first estimate of the spectrum by 
sampling the flux in the three energy bands (0.1-0.5 GeV, 0.5-1 GeV, 1-3 GeV; see Figure 2) 
where the source is clearly detected. A power law fit yields a photon index a — — 1.67 ± 0.25. 
The AGILE spectral points are a factor ~2 below the previous EGRET upper limits (|25l) and 
well above the extrapolation of HESS vF v spectrum to lower energies. The PWN gamma ray 
luminosity in the 0.1-10 GeV band, for a distance of ~290 pc (0 126]), is At\ x 10 33 erg s _1 
corresponding to ~10 -3 E rot (where E rot is the spin-down luminosity of the pulsar). Such a 
luminosity is slightly higher than at VHE energies (9.9 x 10 32 erg s _1 ). 

In the frame of leptonic models, the AGILE measurements are not consistent with a sim- 
ple multiwavelength spectral energy distribution involving a single electron population. The 
AGILE spectral points are one order of magnitude above the fluxes expected from the electron 
population simultaneously fitting synchrotron x-ray emission (peaking at ~1 keV) and inverse- 
Compton (IC) TeV emission §M Ii4)> . 

Additional electron populations should be invoked to explain the observed GeV fluxes. This 
is not surprising in view of the complex morphology of the PWN seen in radio and x-rays, 
where different sites and features of non-thermal emission are present: the anisotropic pulsar 
wind and non-homogeneous SNR reverse shock pressure produce different particle populations 
within the shocked wind. In particular, assuming the same magnetic field (5/iG) reproducing the 
TeV spectral break, the radio synchrotron emitting electrons observed in the Vela X structure 
(1271) may be responsible for the inverse Compton (IC) bump in the GeV band arising from scat- 
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tering on CMBR and Galactic and starlight photon fields, as noticeably predicted by de Jager et 
al. (1771 [7?1) . In fact, the position where AGILE sees the maximum brigthness (Ra. = 08 h 35 m , 
Decl. = — 45°44') is also roughly where the 8.4 GHz radio emission is brightest (|2gl) . AGILE 
data are compatible with the IC parameters modelled by de Jager et al. rfl~5l) (electron spectral 
index 1 .78 and maximum energy ~20 GeV), although our measurements could suggest a higher 
contribution from IC photon seeds. In particular, assuming a starlight energy density of 1.4 eV 
cm -3 and a mean temperature of ~2300 K (1291) . we obtain a good description of the AGILE 
data (Figure 2). 

The AGILE measurements would be incompatible with the scenario of nucleonic gamma- 
ray production in the Vela TeV nebula in the frame of a single primary electron population. 
These models predict very faint GeV emission (< 10 30 erg s _1 ) even when including synchrotron 
and inverse Compton emission from primary and secondary electrons produced by the inelastic 
nuclear scattering (17oT) . On the other hand, the proposed additional electron component scenario 
described above leaves room for uncorrelated GeV- TeV emission, although the comprehensive 
multiwavelength two-component leptonic model (providing strong IC emission on a relatively 
dense photon field) seems to disfavour dominant nucleonic gamma-ray production. In fact, it 
has been found that the thermal particle density at the head of the cocoon, where bright VHE 
gamma-ray emission was found, is a factor of 6 lower than that required by hadronic models 
(UJ). 

The radio emitting region mentioned above appears to be larger (~2x3 square degrees) 
than the AGILE nebula, possibly indicating that IC cooling in the GeV domain is important, 
but we notice that the actual physical size of the GeV nebula could be larger than that presently 
resolved with the available photon statistics, due to the strong Galactic gamma-ray emission 
affecting MeV-GeV energy bands. Instead, the AGILE nebula is similar in shape to the HESS 
nebula, possibly suggesting that the core of HE and VHE emission is produced in the same 
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projected region of Vela X, even if different electron populations are involved. Indeed, different 
spots of bright radio emission (l28l) . possibly associated to electrons injected at different stages 
of pulsar evolution, are embedded within the poorly resolved HE and VHE emission regions. 

High-energy emitting PWNe are thought to be a common phenomenon associated to young 
and energetic pulsars (TJOl) because their IC emission arises mostly from scattering on CMBR 
and starlight fields with no special environmental requirements. On the other hand, PWNe are 
expected to be much weaker than pulsed emission from the associated neutron star, especially 
in the GeV domain where most of the pulsar's spin-down energy is funnelled. 

Indeed, in spite of a PWN gamma ray yield of L^ WN ~ 10~ 3 x E rot , to be compared with the 
typical gamma-ray pulsed luminosity of LP ulsed ~ (10~ 2 -0A)xE ro t, our AGILE observation 
shows that 10 ky old PWNe can match the sensitivities of current GeV instruments. 

Because the gamma-ray luminosity of the PWNe is only a small fraction of the beamed 
emission from the neutron star, the PWN component is difficult to identify in weaker gamma-ray 
pulsars, although it could account for a substantial part of the observed off-pulse flux. However, 
if the beamed emission does not intersect the line of sight to the observer, the PWN component, 
unhindered by the stronger pulsed emission, could be detectable. Energetic pulsars (e.g. E rot ~ 
10 37 erg s _1 ) can power PWNe with gamma-ray luminosities matching the flux (~10~ 8 -10~ 7 
ph cm" 2 s -1 ; E > 100 MeV) of a class of unidentified EGRET sources (|22|) . as well as a subset 
of the ones detected by AGILE and Fermi (TJ/I |32l) . when placed within few kiloparsecs. The 
roughly isotropic emission from such undisturbed PWNe would not yield pulsations and, as 
a class, they could contribute to the population of Galactic unidentified sources still awaiting 
multiwavelength association (TJillTJl) . 
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Figure 1: Top: Gaussian-smoothed AGILE intensity map (ph cm~ 2 s _1 sr _1 with 0.25x0.25 
square degrees pixel size) at E > 400 MeV around the Vela pulsar, including only off-pulse 
events (i.e. discarding events with pulsar phase corresponding to Vela pulsed emission); the neu- 
tron star position is marked with the black cross, green circles are the 68% confidence contours 
for the position of AGL J0848-4242 and AGL J0834-4539 [Vela X]. The AGILE E > 400 MeV 
energy band is well suited for gamma-ray imaging and provides a good compromise between 
the instrument effective area (~400 cm 2 ; ~100 counts from AGL J0834-4539) and point spread 
function (~1°, 68% containment radius), both parameters decreasing with energy. Bottom: the 
gamma-ray diffuse source AGL J0834^4539. AGILE contours (bottom — left) are in the range 
(1.4-1.6) x 10~ 4 ph cm -2 s _1 sr -1 , with step 4xl0 -6 . HESS contours (bottom — right) arc 
taken from (17D1) . 
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Figure 2: Gamma-ray high- and very-high energy spectral distribution (yF v ) of the Vela X 
PWN. HESS data fit an IC process (scattering on CMBR) related to electron power law index 
2.0 with a break at 67 TeV and a total energy content of 2.2 x 10 45 erg ((721). AGILE data 
are compatible with IC emission from the additional electron component well reproducing the 
observed total radio spectrum (E tot = 4 x 10 48 erg) assuming the same ~ 5 /xG field strength 
as required by the TeV spectral break. Unlike the TeV IC emission, GeV IC scattering is 
within Thomson limit. Thus, in addition to the CMBR component (photon density n p h = 0.25 
eV cm" 3 , photon energy E ph = 1(T 3 eV), also FIR (n ph = 0.3 eV cm" 3 , E ph = 1(T 2 eV) 
and starlight (n p h = 1.4 eV cm -3 , E ph = 1 eV) photon fields can significantly contribute to 
the high-energy IC counterpart of the radio spectrum fitting AGILE data (dot — dashed lines: 
CMBR (a), FIR (6), starlight (c); thick line: total IC spectrum). 
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